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ABSTRACT 

We report the discovery of a sub-Jupiter mass exoplanet transiting a magnitude V = 11.7 host star 1SWASP J030928. 54+304024.7. 
A simultaneous fit to the transit photometry and radial-velocity measurements yield a planet mass M v = 0.53 + 0.07 Mj, radius 



no 1+0.061 

yl -0.03 ' 



p - xj.si _ 003 Rj and an orbital period of 3. 722465 ^omooos days- The host star is of spectral type K3V, with a spectral analysis yielding 
an effective temperature of 4800 + 100 K and logg = 4.45 ± 0.2. It is amongst the smallest, least massive and lowest luminosity stars 
known to harbour a transiting exoplanet. WASP-llb is the third least strongly irradiated transiting exoplanet discovered to date, 
experiencing an incident flux F p = 1.9 x 10 s ergs -1 cirT 2 and having an equilibrium temperature T cq \ = 960 ± 70 K. 



1. Introduction 

Observations of planets that transit their host star represent 
the current best opportunity to test models of the internal 
.£h ! structure of exoplanets and of their formation and evolution. 
. Since the first detection of an exoplanetary transit signature 
Uh ' dCharbonneau et all l2Q0Ct iHenrv et aTJ |200Q|) over fifty tran- 
siting planetary systems have been identified. A number of 
wide-field surveys are in progress with t he goal of detectin g 
trans iting exoplanets, for exam ple OGLE ( Udalski et al. 2002), 
XO dMcCullough et alJ |2005h. HAT dBakos et alj|2004l). TrES 
dO'Donovan et al.ll2006l) and WASP dPollacco et alj|2006l) . 

The WASP project operates two identical instruments, at La 
Palma in the Northern hemisphere, and at Sutherland in South 
Africa in the Southern hemisphere. Each telescope has a field of 
view of just under 500 square degrees. The WASP survey is sen- 
sitive to planetary transit signatures in the light-curves of hosts 
in the magnitude range V~9-13. A detailed description of the 
telescope har dware, observing stra tegy and pipeline data analy- 
sis is given in lPollacco et alj (|2006). 

In this paper we report the discovery of WASP-llb, a 
sub-Jupiter mass gas giant planet in orbit about the host star 
1SWASP J030928.54+304024.7. We present the WASP discov- 
ery photometry plus higher precision optical follow-up and ra- 
dial velocity measurements which taken together confirm the 
planetary nature of WASP- 1 lb. 



2. Observations 

2.1. WASP photometry 

The host star 1SWASP J030928.54+304024.7 (= USNO- 
B1.0 1206-0003989 = 2MASS 03092855+3040249; hereafter 
labelled WASP-11) was observed by SuperWASP-N during the 
2004, 2006 and 2007 observing seasons, covering the inter- 
vals 2004 July 08 to 2004 September 29, 2006 September 
09 to 2007 January 20 and 2007 September 04 to 2007 
December 12 respectively. The pipeline-processed data were de- 
tr ended and searched for tran sits using the methods described 
in IColher Cameron et alJ d2006l) . yielding a detection of a peri- 
odic transit-like signature with a period of 3.722 days. A total of 
ten transits are observed in data from all three observing seasons 
(TableQ] FigureQ). 



2.2. Photometric follow-up 

WASP- 1 1 was followed-up with the 2-m Liverpool telescope on 
La Palma as part of the Canarian Observatories' International 
Time Programme for 2007-08. We used the 2048 x 2048 pixel 
EEV CCD42-40 imaging camera giving a scale of 0.27 arcsec- 
onds/pixel in 2x2 bin mode and a total field of view of ~ 4.6x4.6 
arcminutes 2 . Observations were taken during the transit of 2008 
January 14, and consist of 656 images of 10 seconds exposure 
in the Sloan z band. The night was non-photometric and with 
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Table 1. WASP-N survey coverage of WASP- 1 1 



Table 2. Radial velocity measurements of WASP- 1 1 



Season 


Camera 


N pts 


N tr 


To 

BJD-2400000.0 


P 

(days) 


2004 


103 


1756 


4 


53240.921696 


3.7220 


2006 


144 


2679 


3 


54056.140758 


3.7223 


2007 


146 


2750 


2 


54346.4883 


3.7226 


2007 


147 


729 


1 








Fig. 1. SuperWASP-N photometry of WASP- 1 1 from the 2004, 
2006 and 2007 seasons. The data have been de-trended using 
the sysrem scheme described in (ICollier Cameron et alj [2006) 
and are plotted here phase-folded on the best-fit period from the 
MCMC analysis (section|3]l. 



seeing varying from 0.9 to 2.2 arcsec during the four hour long 
observing run. 

The images were bias subtracted and flat-field corrected with 
a stacked twilight flat-field image. All the science images were 
also corrected for the fringing effect. The autoguider did not 
work during our observations and a maximum positional shift 
of 17.5 arcsec of the stars within the frame was noticed. After 
aligning the images with respect to the first target image, aper- 
ture photometry were performed around the target and compari- 
son stars using an aperture of 20 pixels (5'.'4) radius. Three bright 
non-variable comparison stars were available in the target field 
with which to perform differential photometry. 

Further observations of WASP- 11 were made with the 
Keele University Observatory 60 cm Thornton Reflector on 2008 
February 09 and 13. This telescope is equipped with a 765 x 510 
pixel Santa Barbara Instrument Group (SBIG) ST7 CCD at the 
f/4.5 Newtonian focus, giving a 0.68 arcsecond/pixel resolution 
and a 8.63 x 5.75 arcminute field of view. Conditions were pho- 
tometric throughout both nights, although the transit of February 
9 ended at an airmass of 4 and cryogenics problems on the night 
of February 13 may have led to some frosting on the CCD dewar 
window during the first few exposures. Tracking errors and spu- 
rious electronic noise mean that systematic noise is introduced 
into the system at an estimated level of 4 millimag with period- 
icities of 2 (worm error) and 20 minutes (presently of unknown 
origin). No corrections have been applied for these effects. 

Altogether (237 + 276)x 30 s observations in the R band 
were obtained. After applying corrections for bias, dark cur- 
rent and flat fielding in the usual way, aperture photometry on 
WASP- 11 and the comparison star USNO-B1.0 1207-0040657 
were performed using the commercial software AIP4Win 
dBerrv & Burnelll 12005). The resulting lightcurves from both 
Liverpool Telescope and Keele 60 cm observations (Figure|2]top 
panel) confirm the presence of a transit. 
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0.000 


SOPHIE 
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SOPHIE 
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0.0076 


-0.002 


SOPHIE 
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0.0106 


-0.002 


SOPHIE 
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0.0143 


-0.035 


SOPHIE 


2454512.3848 


4.9482 


0.0096 


0.022 


SOPHIE 



2.3. Radial velocity follow-up 

Initial spectroscopic observations were obtained using the Fibre- 
fed Echelle Spectrograph (FIES) mounted on the 2.5-m Nordic 
Optical Telescope. A total of five radial velocity points were 
obtained during 2007 December 27-31 and 2008 January 25. 
WASP- 1 1 was observed with an exposure time of 1800s giving a 
signal-to-noise ratio of around 70-80 at 5500A. FIES was used 
in medium resolution mode with R=46000 with simultaneous 
ThAr calibration. We used the bespoke data reduction package 
FIEStool to extract the spectra and a specially developed IDL 
line-fitting code to obtain radial velocities with a precision of 
20-25ms- 1 . 

Radial velocity measurements of WASP- 1 1 were also made 
with the Observatoire de H aute-Provence's 1.93 m telescope an d 
the SOPHIE spectrograph dBouchv & The Sophie Te am 2006), 
over the 8 nights 2008 February 11 - 15; a total of 7 usable 
spectra were acquired. SOPHIE is an environmentally stabilized 
spectrograph designed to give long-term stability at the level of 
a few ms _1 . We used the instrument in its medium resolution 
mode, acquiring simultaneous star and sky spectra through sep- 
arate fibres with a resolution of R=48000. Thorium- Argon cal- 
ibration images were taken at the start and end of each night, 
and at 2- to 3-hourly intervals throughout the night. The radial- 
velocity drift never exceeded 2-3 m s , even on a night-to-night 
basis. 

Conditions during the SOPHIE observing run were photo- 
metric throughout, though all nights were affected by strong 
moonlight. Integrations of 1080 s yielded a peak signal-to- 
noise per resolution element of around ~30^40. The spectra 
were cross-correlated against a K5V template provided by the 
SOPHIE control and reduction software. 

In all SOPHIE spectra the cross-correlation functions (CCF) 
were contaminated by the strong moonlight. We corrected them 
by using the CCF from the background light's spectrum (mostly 
the Moon) in the sky fibre. We then scaled both CCFs using the 
difference of efficiency between the two fibres. Finally we sub- 
tracted the corresponding CCF of the background light from the 
star fibre, and fitted the resulting function by a Gaussian. The pa- 
rameters obtained allow us to compute the photon-noise uncer- 
tainty of the corrected radial velocity measurement (cr RV ), using 
the relation 

cr RV = 3.4 ^(FWHM)/(S/N x Contrast) 

Overall our SOPHIE RV measurements have an average photon- 
noise uncertainty of 10.3 m s _1 . The measured barycentric radial 
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Fig. 2. The best-fit model from the simultaneous MCMC fit to 
the available photometry (top panel) and radial velocity data 
(lower panel). The fitted zero-point offset between the NOT and 
SOPHIE radial-velocity measurements (5.4 + 0.4m s -1 ) is re- 
moved in this plot. 



velocity (Table|2] Figure[2]lower panel) show a sinusoidal varia- 
tion of half-amplitude ~ 90 ms -1 about a centre-of-mass RV of 
~ 4.9kms _1 , consistent with the presence a companion of plan- 
etary mass. The period and ephemeris of the RV variation are 
consistent with those of found by the transit search. 

An analysis of the line-bisector spans shows no significant 
correlation with radial velocity (Figure[3j, as would be expected 
if the observed radial velocity variations w ere due to a dilute d 
eclipsing binary or chromospheric activity (Oue loz et al.l l2001). 




4.9 4.95 
Radial velocity (km/s) 



Fig. 3. The line bi-sector against velocity for WASP-1 1, showing 
no evidence of correlation. 



Table 3. System parameters of WASP-1 1 derived from a simul- 
taneous MCMC analysis of the available photometric and radial- 
velocity measurements. Quoted uncertainties define the lcr con- 
fidence intervals. 



Transit epoch (HJD), To 


2454473.05586 ±0.0002 




Orbital period, P 


■2 7994£c+0.000006 
— 0iX)0008 


days 




0162+ 0003 
u.uiuz. _ 00002 




Transit duration 


2 556 +a029 


hours 


Impact parameter, b 


054 +016S 


R, 


Reflex velocity, K] 


0.0821 ± 0.0074 


kms~' 


Centre-of-mass velocity, y 


4.9077 + 0.0015 


kms" 1 


Orbital eccentricity, e 


= 0.0 




Orbital inclination, (' 


89.8*3 


deg 


Orbital separation, a 


0.043 + 0.002 


AU 


Stellar mass, M t 


77+ 010 

' ' -0.08 




Stellar radius, 


74+ - 04 
-Q.Q3 


Ro 


Planet radius, R p 


91+0-05 


R, 


Planet mass, M p 


0.53 ± 0.07 


Mj 


l°ggp (cgs) 


3 16 +a04 




Planet density, p p 


69 +0 07 


Pj 


Planet T cql (A=0;f=l) 


960 ± 70 


K 



3. System parameters 

3.1. Stellar parameters 

In order to perform a detailed spectroscopic analysis of the stel- 
lar atmospheric properties of WASP-1 1, we merged the avail- 
able FIES spectra into one high-quality spectrum, carefully re- 
moving any radial velocity signature during the process. This 
merged spectrum was then continuum-normalized with a very 
low order polynomial to retain the shape of the broadest spectral 
features. The total signal-to-noise of the combined spectrum was 
around 200 per resolution element. We were not able to include 
the SOPHIE spectra in this analysis, because these spectra were 
obtained with the HE (high-efficiency) mode which is known to 
suffer from problems with removal of the blaze function. 

For our analysis we followed the same pr ocedure as for 
the sp ectroscopic ch aracterization of WA SP-1 (Stempels et al.l 
120071) and WASP-3 dPollacco et a lJ EooH). We used the pack - 
age Spectroscopy Made Easy (SME. lValenti & Piskunov| [l996). 
which combines spectral synthesis with multidimensional x 1 
minimization to determine which atmospheric parameters best 
reproduce the observed spectrum of WASP-1 1 (effective tem- 
perature T e ff, surface gravity \ogg, metallicity [M/H], projected 
radial velocity v sin i, systemic radial velocity v rad , microturbu- 
lence v m ; c and the macro turbulence v mac ). For a more detailed 
descr iption of the spectral s ynthesis and our assumptions we re- 
fer to lStempels et al.l d2007l) . 

The four spectral regions we used in our analysis are (1) 
5160-5190A, covering the gravity-sensitive Mg b triplet (2) 
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Wavelength (A) 

Fig. 4. A comparison between the observed FIES spectrum of WASP- 1 1 and the calculated spectrum obtained from spectral syn- 
thesis with SME. The white regions are excluded from the spectral analysis, mainly because of the presence of telluric absorption. 
Light shaded regions were used to determine the continuum level, and the remaining dark shaded regions to determine the stellar 
atmospheric parameters. 



5850-5950A, with the temperature and gravity-sensitive Na 
i D doublet (3) 6000-621 OA, containing a wealth of differ- 
ent metal lines, providing leverage on the metallicity, and (4) 
6520-6600A, covering the strongly temperature-sensitive H- 
alpha line. A comparison between the observed FIES spectrum 
and the synthetic spectrum is shown in Figure [4] The spectral 
analysis yields an effective temperature T e ff = 4800 + 100K, 
logg = 4.45 + 0.2, [M/H] = 0.0 + 0.2 and vsini < 6.0km s -1 . 
These parameters correspond to spectral type of K3V. A close 
examination of the region around the Li i 6708 shows no evi- 
dence of such a feature, suggesting that the lithium abundance is 
very low. 

3.2. Planet parameters 

To determine the planetary and orbital parameters the SOPHIE 
and NOT FIES radial velocity measurements were combined 
with the photometry from WASP and the Liverpool Telescope 
in a simultaneous fit using the Markov Chain Monte Carlo 
( MCMC) technique. T he details of this process are described 
in lPollacco et al.l 12008). An initial fit showed that the orbital ec- 
centricity (e - 0.086*Q Qg2) was poorly constrained by the avail- 
able data and nearly consistent with zero. We therefore fixed the 
eccentricity parameter at zero in a further fits. Figure [2] shows 



the best-fitting models. The best-fit parameters (Table [3]) show 
WASP- 1 lb to have a mass M = 0.53 + 0.07 Mj and a radius of 
R = 0.91!<J$ Rj. 



4. Discussion 

The system parameters derived here place WASP- lib towards 
the lower end of the mass range of known transiting planets, 
falling approximately mid-way between the masses of Jupiter 
and Saturn. The host star WASP-11 is also amongst the small- 
est and lowest luminosity stars known to host a transiting 
planet, however it is relatively nearby and thus quite bright 
(V = 11.7). WASP- lib is irradiated by a stellar flux F p = 
1.9 x 10 8 ergcrrr 2 s _1 at the sub-stellar point making it the 
third least heavily irradiated transiting planet after GJ436b and 
HD 17 156b. We compute an equilibrium temperature for WASP- 
11b of r eq i(A = 0;/ = 1) = 960 + 70 K, which makes it more 
typical of the bulk of known exoplanets than of the "hot Jupiter" 
class most commonly found by the transit method. 

Theoretical models of the atmospheres of hot giant exoplan- 
ets dFortnev et al.l 120061: iBurrows etai1l2007l) have shown that 
heavy irradiation can lead to the development of a temperature 
inversion and a hot stratosphere. This is due to the absorption 
of stellar flux by an atmospheric absorber, possibly TiO and 
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Fig. 5. The position of WASP-1 1 in the R/M 1/3 - T eS plane 
Evolu tionary tracks for a solar metallici t y sta r from Baraffe ^t al.l 
( 1998) (upper panel) and lGirardi et al.l (l2000h (lower panel) are 
plotted along with isochrones for ages lOMyr (solid), 1 Gyr 
(dashed), 5 Gyr (dot-dashed), 10 Gyr (dotted). Evolutionary 
mass tracks are shown for 0.7, 0.8, 0.9 and 1.0 M . 



VO. In both sets of models the magnitude of the incident stel- 
lar flux is the key controlling variable determining whether a 
given extra-solar giant planet (E GP) will posse s s a h ot strato- 
sphere. Recent observations by Mac halek et al.l d2008l) of sec- 
ondary transits of XO-lb using the Spitzer Space Telescope 
suggest the presence of a temperature inversion in the atmo- 
sphere of that exoplane t. On the other hand ana logous obser- 
vations of HD 189733b dCharbonneau et alj|2008h show no evi- 
dence for an inversion, despite the irradiating fluxes of XO-lb 
and HD189733b being almost identical (F p = 0.49 x 10 9 and 
F p = 0.47 x 10 y erg cirT 2 s -1 respectively). This strongly sug- 
gests that the incident stellar flux is not the sole controlling pa- 
rameter determining the presence of the inversion, a likelihood 
which the authors of the atmosphere models readily point out 
themselves. Further observations of planets particularly in the 
low-irradiation regime are required to help parameterise the ther- 
mal inversion. WASP-1 lb is amongst the nearest and brightest 
low-irradiation EGPs making it a good candidate for such stud- 
ies. Moreover we note that the orbital eccentricity of WASP- lib 
is much lower than the other two bright low-irradiation transit- 
ing exoplanets, GJ436b and HD 17 156b (e = 0.15 and e = 0.67 
respectively). As a consequence the secular variation in irradi- 
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Fig. 6. Planetary mass-radius relations as a functio n of core mass 
and sy stem age, interpolated from the models of iFortnev et al.l 
J2007T) . 



ation around the orbit will be correspondingly lower in WASP- 
11b, removing a potentially complicating factor when compar- 
ing follow-up observations with predictions from atmospheric 
models developed assuming steady-state irradiation. 

To estimate the age of the WASP-1 1 we compared the ob- 
served stellar density and temperature against t he evolution- 
ary mo dels of low- and interm ediate-mass stars of iGirardi et al.l 
(120001) and lBaraffe et al l d!998l) . In Figure|5]we plot the position 
of WASP-1 1 in the R/M 1 ^ versus r e ff plane atop isochrones of 
different ages from the two models. For such a cool star, the 
isochrones are closely spaced in this parameter plane due to the 
slow post-main-sequence evolution of late-type stars. The sets of 
isochrones from the two models overlap in this regime, and both 
models suggest the same mass and age for the host star. WASP- 
1 1 falls above the 10 Gyr isochrone for both models, though it is 
consistent with this age within the errors. The very low lithium 
abundance also points to ward WASP-1 1 being ^ 1-2 Gyr old 
(ISestito & Randichll2005l) . We inv estigated using gyrochronol- 
ogy to age the host star, following iBarnesI (120071) . however we 
were unable to measure a definite rotational period. No rotation 
modulation was detected in the lightcurve to an amplitude limit 
of a few milli-magnitudes. The spectral analysis furnishes only 
an upper-limit to v sin i, so no rotational period can be deter- 
mined in that way. Taken together these factors are all consis- 
tent with WASP-1 1 being an old star, older than maybe 1 Gyr, 
however it is not possible to be more definite than that with the 
availab le data. 

IFortnev et al.l d2007) present models of the evolution of plan- 
etary radius over a range of planetary masses and orbital dis- 
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tances, and under the assumption of the presence of a dense core 
of various masses up to 100 M®. To compare our results with 
the Fortney et al. models we plotted the modelled mass-radius 
relation as a function of core mass in Figure [6] To account for 
the lower-than-Solar luminosity of the host star WASP- 11 we 
calculated the orbital distance a Q = a(M*/M ) -3 - 5 ' 2 at which a 
planet in orbit about the Sun would receive the same incident 
stellar flux as WASP- lib does from its host. We then interpo- 
lated the models of Fortney et al. to this effective orbital distance 
(a = 0.068 for WASP-llb). As the age of the WASP- 11 sys- 
tem is poorly constrained we compare our results with the mod- 
elled mass-radius relation at 300 Myr, 1 Gyr and 4.5 Gyr. We find 
that the radius of WASP- 1 lb is consistent with the presence of a 
dense core with a mass in the range M core ~ 42-77 M e for a sys- 
tem age of 300 Myr, M core -33-67 M® at 1 Gyr, and M core -22- 
56 M® at 4.5 Gyr. 

Acknowledgements. The WASP Consortium consists of astronomers primarily 
from the Queen's University Belfast, Keele, Leicester, The Open University, and 
St Andrews, the Isaac Newton Group (La Palma), the Instituto de Astrofisica 
de Canarias (Tenerife) and the South African Astronomical Observatory. The 
SuperWASP-N and WASP-S Cameras were constructed and operated with 
funds made available from Consortium Universities and the UK's Science 
and Technology Facilities Council. SOPHIE observations have been funded 
by the Optical Infrared Coordination network (OPTICON), a major interna- 
tional collaboration supported by the Research Infrastructures Programme of the 
European Commission's Sixth Framework Programme. FIES observations were 
made with the Nordic Optical Telescope, operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de 
Canarias. We extend our thanks to the Director and staff of the Isaac Newton 
Group of Telescopes for their support of SuperWASP-N operations, and the 
Director and staff of the Observatoire de Haute-Provence for their support of 
the SOPHIE spectrograph. The Liverpool Telescope is operated on the island of 
La Palma by Liverpool John Moores University in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias with finan- 
cial support from the UK Science and Technology Facilities Council. 



References 

Bakos, G., Noyes, R. W., Kovacs, G, et al. 2004, PASP, 1 16, 266 

Baraffe, I., Chabrier, G., Allard, F, & Hauschildt, P. H. 1998, A&A, 337, 403 

Barnes, S. A. 2007, ApJ, 669, 1167 

Berry, R. & Burnell, J. 2005, The handbook of astronomical image processing 
(The handbook of astronomical image processing, 2nd ed., by R. Berry and 
J. Burnell. xxviii, 684 p., 1 CD-ROM (incl. Astronomical Image Processing 
Software AIP4WIN, v.2.0). Richmond, VA: Willmann-Bell, 2005) 

Bouchy, F. & The Sophie Team. 2006, in Tenth Anniversary of 51 Peg-b: 
Status of and prospects for hot Jupiter studies, ed. L. Arnold, F. Bouchy, & 
C.Moutou, 319-325 

Burrows, A., Hubeny, I., Budaj, J., Knutson, H. A., & Charbonneau, D. 2007, 
ApJ, 668, L171 

Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, M. 2000, ApJ, 529, 
L45 

Charbonneau, D., Knutson, H. A., Barman, T, et al. 2008, ArXiv e-prints, 802 
Collier Cameron, A., Pollacco, D., Street, R. A., et al. 2006, MNRAS, 373, 799 
Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 
Fortney, J. J., Saumon, D., Marley, M. S., Lodders, K, & Freedman, R. S. 2006, 
ApJ, 642, 495 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 
Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, ApJ, 529, L41 
Machalek, P., McCullough, P. R., Burke, C. J., et al. 2008, ArXiv e-prints, 805 
McCullough, P. R., Stys, J. E., Valenti, J. A., et al. 2005, PASP, 117, 783 
O'Donovan, F. T., Charbonneau, D., Mandushev, G., et al. 2006, ApJ, 651, L61 
Pollacco, D., Skillen, I., Collier Cameron, A., et al. 2008, MNRAS, 385, 1576 
Pollacco, D. L., Skillen, I., Cameron, A. C, et al. 2006, PASP, 118, 1407 
Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279 
Sestito, P. & Randich, S. 2005, A&A, 442, 615 

Stempels, H. C, Collier Cameron, A., Hebb, L., Smalley, B., & Frandsen, S. 

2007, MNRAS, 379, 773 
Udalski, A., Paczynski, B., Zebrun, K, et al. 2002, Acta Astronomica, 52, 1 
Valenti, J. A. & Piskunov, N. 1996, A&AS, 118, 595 



